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ABSTRACT 

PSR J1802— 2124 is a 12.6-ms pulsar in a 16.8-hour binary orbit with a relatively massive white 
dwarf (WD) companion. These properties make it a member of the intermediate-mass class of binary 
pulsar (IMBP) systems. We have been timing this pulsar since its discovery in 2002. Concentrated 
observations at the Green Bank Telescope, augmented with data from the Parkes and Nangay observa- 
tories, have allowed us to determine the general relativistic Shapiro delay. This has yielded pulsar and 
white dwarf mass measurements of 1.24 ± 0.11 M© and 0.78 ± 0.04 M© (68% confidence), respectively. 
The low mass of the pulsar, the high mass of the WD companion, the short orbital period, and the 
pulsar spin period may be explained by the system having gone through a common-envelope phase in 
its evolution. We argue that selection effects may contribute to the relatively small number of known 
IMBPs. 

Subject headings: pulsars: general pulsars: individual (PSR J1802— 2124) binaries: general stars: 
evolution 



1. INTRODUCTION 

In the standard theory of pulsar spin-up, a neutron 
star (NS) in a binary system accretes matter from its 
companion star. This serves to transfer angular momen- 
tum to the NS, increasing the spin frequency of the pul- 
sar (e.g., Alpar et al. 1982). The type and duration of 
mass transfer onto the pulsar determines the final spin 
period and depends a great deal on the nature and evo- 
lution of the system (for reviews of binary pulsar systems 
and their evolution; e.g., Bhattacharya & van den Hcuvcl 
1991; Phinney & Kulkarni 1994; Stairs 2004; Tauris & 
van den Heuvel 2006). 
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The measured mass distribution of NSs in pulsar bi- 
nary systems is more diverse than previously thought. 
Observations show that many pulsars have masses which 
lie significantly outside the oft-cited statistical average 
of 1.35 ±0.04 M G (Thorsett & Chakrabarty 1999). This 
highlights the need to invoke a variety of evolutionary 
scenarios in order to explain the collection of observed 
pulsar binary systems. Most double-neutron-star (DNS) 
binary systems are thought to have undergone common- 
envelope (CE) evolution, in which the neutron star re- 
sulting from the first supernova spirals into the envelope 
of the companion star; the CE is subsequently expelled 
from the system (e.g., Tauris & van den Heuvel 2006). 
That this mass-transfer stage in the evolution of these 
systems is short-lived is evidenced in part by the rela- 
tively long ~ 10 — 100-ms pulsar spin periods observed 
in DNS binaries. 

The majority of known neutron star-white dwarf (NS- 
WD) binaries have pulsars that spin with millisecond pe- 
riods (< 10 ms) and have very low eccentricities. Such 
short rotation periods indicate that the pulsars have un- 
dergone a relatively long, stable period of accretion of 
material from the companions' outer envelopes, during 
which they are seen as low-mass X-ray binaries (Bhat- 
tacharya & van den Heuvel 1991; Wijnands & van der 
Klis 1998; Archibald et al. 2009). In the process, the 
matter-donating stars lose an appreciable amount of 
mass. This conclusion is supported by the relatively low 
masses of the WDs found in these binaries, referred to as 
low-mass binary pulsar (LMBP) systems (e.g., Tauris & 
van den Heuvel 2006). 

In contrast to the LMBP binaries are the intermediate- 
mass binary pulsar class (IMBPs; Camilo et al. 1996; Ed- 
wards & Bailes 2001a). These systems are characterized 
by pulsar spin periods of tens of milliseconds, massive 
carbon-oxygen (CO) or oxygen-neon-magnesium (O-Nc- 
Mg) WD companions (> 0.4 M©), and orbital eccentric- 
ities which, while still small, are significantly larger than 
those of LMBP systems. 
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TABLE 1 

Summary of observations of PSR J1802— 2124 



Telescope 


Instrument 


Center 


Total effective 


Integration 


Number 


Date span 






frequency 


bandwidth 


time 


of TOAs 


(MJD) 






(MHz) 


(MHz) 


(min) 






Parkes 


Filtcrbank 


1390 


256 


typically 20 


106 


52605-54910 


GBT 


GASP 


1400 


64-96 


3 


2233 


53441-54950 


Nangay 


BON 


1398 


64 - 128 


14 


107 


54188-54806 



Several formation scenarios have been suggested to ex- 
plain the existence of IMBPs (e.g., Li 2002). One idea 
is that, as with DNS binaries, the NS spirals into the 
envelope of its companion to form a CE which is then 
promptly ejected from the system (e.g., van den Heuvel 
1994). This is supported by the short orbital periods (Pb) 
seen in many of the IMBP systems. It has also been pro- 
posed, however, that a neutron star within the envelope 
of its companion will be forced to undergo hypercritical 
accretion, becoming a black hole and thus rendering the 
system unobservable (e.g., Chevalier 1993; Brown 1995; 
Brown et al. 2001, hereafter B01). 

Tauris et al. (2000) argued that systems with heavy 
CO WD companions and orbits with 3 < Pb < 70 days 
can undergo, and survive, a short-lived phase of highly 
super-Eddington mass transfer to the NS. Here, the in- 
spiral that results in a CE is avoided if the re-radiated 
accretion energy is great enough to evaporate most of 
the transferred material before it approaches the NS too 
closely (see also Taam et al. 2000). Still, this scheme 
docs not work for IMBPs with Pb < 3 days, suggesting 
the need to invoke CE evolution to explain their exis- 
tence. B01 put forward a possible alternate formation 
scenario in which the two progenitors are main sequence 
stars of similar mass, which evolve to form two helium 
cores. This is similar to a related scenario for the evolu- 
tion of close DNS binaries (Brown 1995), the difference 
being that to form an IMBP, one of the stars would be 
just below the mass threshold for NS formation, becom- 
ing a WD instead. It is clear that the evolution of IMBPs 
remains an open question. 

There are now sixteen known IMBP systems. Until 
now, only one of these, PSR J0621+1002, has a pulsar 
with a measured mass (mi = 1.70±8;^ M ; Nice et al. 
2008) , which is somewhat larger than the range of masses 
seen in well-measured DNS systems (e.g., Stairs 2008). 
Here we present results from timing of PSR J1802— 2124, 
another member of the class of IMBPs, discovered in the 
Parkes Multibcam Pulsar Survey (Faulkner et al. 2004). 
In the case of PSR J1802— 2124, its heavy companion led 
us to believe that the system would be a good candidate 
for measuring the Shapiro delay of the pulses in the grav- 
itational potential of the white dwarf. We were, in fact, 
able to measure such an effect, enabling us to determine 
the individual masses of each member of the binary sys- 
tem. In this article, we discuss these measurements and 
their implications for reconstructing the formation and 
evolutionary histories of this system and others like it 1 . 

2. OBSERVATIONS 

1 A preliminary version of these results was reported in Fcrdman 
(2008) 



We used three observatories to collect pulsar data. In 
what follows, we describe the instruments used and the 
data obtained at each of these telescopes. A summary of 
the observing details is found in Table 1. 

2.1. Parkes 

We have added to the data set presented by Faulkner 
et al. (2004) using the 64-m Parkes telescope in Australia. 
Observations were carried out at regular intervals using a 
2 x 512 x 0.5-MHz filterbank centered at 1390 MHz, each 
typically 20 minutes in duration. The data from each 
channel were detected and the two polarizations summed 
in hardware before 1-bit digitization every 80 — 250 [is. 
The data were recorded to tape and subsequently folded 
off-line. Parkes data used for this work were collected 
at 106 epochs over 6.3 years. This long timing baseline 
was particularly useful for measurements of astrometric 
parameters. Discovery and some initial timing data were 
taken with a different filterbank (3-MHz channels) and 
were excluded from our analysis. 

2.2. Green Bank 

We also used the 100-m Robert C. Byrd Green Bank 
Telescope (GBT) in West Virginia. Data-taking at the 
GBT was performed with the Green Bank Astronomical 
Signal Processor (GASP; Demorest 2007). GASP is a 
flexible baseband system, which performs 8-bit Nyquist 
sampling of the incoming data stream at 0.25/xs inter- 
vals in both orthogonal polarizations. The signal was 
divided into 16 or 24 x 4-MHz channels 2 centered near 
1400 MHz. The incoming data stream was then coher- 
ently dedispersed (Hankins & Rickett 1975) in software. 
After this, the signals were folded at the pulse period to 
form pulse profiles, typically representing 3-minute in- 
tegration times. These were usually flux-calibrated in 
each polarization by using as a reference the signal from 
a noise diode source that was injected at the receiver. 
When calibration data were not available, we normalized 
the profile data in each polarization by the root-mcan- 
square (rms) value of the corresponding off-pulse signal. 
The data were finally summed over both polarizations 
and across all frequency channels to give the total power 
signal (for further details on GASP operation and data 
reduction, see Demorest 2007; Ferdman 2008). 

Using the GBT, we obtained a total of 42 epochs of 
data that span more than 4 years. Ten of these obser- 
vations consisted of observing sessions between 6 and 8 
hours long. These were scheduled so as to sample as fully 
as possible the orbit of the system, of particular impor- 
tance for the detection of Shapiro delay. The rest of the 

2 The number of channels used occasionally varied due to radio 
frequency interference and available computing resources. 
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Fig. 1— Template profile for PSR J1802-2124. Top: High 
signal-to-noise profile created by all useable data taken at the GBT 
with the GASP pulsar backend instrument. Bottom: Profile cre- 
ated by fitting multiple Gaussians to the PSR J1802— 2124 pulse 
profile. 



data consisted of approximately 1-2 hours of observing. 
These GBT observations comprised the vast majority of 
our data set (65% by time, 97% by weight). 

2.3. Nangay 

We included in our data set observations of 
PSR J1802-2124 taken by the 94-m circular-dish equiv- 
alent Nangay telescope in France. These data were 
recorded with the Berkeley-Orleans-Nangay (BON) pul- 
sar backend, a sister system to the GASP instrument 
at the GBT. The BON instrument is also a baseband 
recorder, which performs coherent dedispersion on the 
incoming data stream in real time. The data originally 
consisted of 16 x 4-MHz channels, and since 2008 July 
25 has been increased to include 32 frequency channels, 
centered at 1398 MHz in both cases. As with the GASP 
backend, the data were detected and folded after sig- 
nal dedispersion was performed. Flux calibration was 
not available for the Nangay data, and so we normal- 
ized each hand of polarization by its off-pulse rms signal 
before obtaining total power pulse profiles. This is ac- 
ceptable since Nangay is a meridian-style telescope, and 
observed PSR J1802— 2124 for approximately 1 hour per 
day; during this time, the telescope gain was not ex- 
pected to change significantly. In all, we collected data 
at 26 epochs over 1.7 years. The output profiles repre- 
sent summation across the observing bandwidth with a 
typical integration time of 14 minutes, or approximately 
4 scans per observing day. Although the Nangay data set 
is small compared with that of the GBT data, it proved 
useful in filling an observational gap between 2007 April 
18 and 2007 October 2007, during which time the GBT 
was undergoing track repair. 

3. TIMING ANALYSIS 

In order to determine the pulse times-of-arrival 
(TOAs), we constructed a template profile for 
PSR J1802-2124. This was done by first aligning in 
phase all individual GBT-derived 1400-MHz pulse pro- 
files, then computing a simple summation of the data, 
weighted by the noise in the off-pulse regions of the in- 
put profiles. Finally, we performed a multiple-Gaussian 



fit to this high signal-to-noise profile, obtaining a zero- 
noise reference template (Kramer et al. 1994; Kramer 
1994), shown in Figure 1. This method allowed us to 
eliminate low-level (< 1% of the peak height) structure 
that was seen to either side of the main pulse in the cu- 
mulative profile; we suspect that this was likely due to 
aliasing that occurs when the channel filters do not cut 
off sharply at the Nyquist frequency (see, e.g., Shrauner 
1997). This noise-free template was then used to calcu- 
late TOAs from both GBT and Nangay data, as their 
hardware configuration and output profile data format 
are similar. The same template was also applied to the 
Parkes data, but was first computed with 256 bins, then 
convolved with a top-hat function representing the in- 
channel dispersive smearing at 1390 MHz, then reduced 
to 64 bins. This was done using the psrchive software 
(Hotan et al. 2004). 

Pulse TOAs were then calculated by cross-correlating 
each pulse profile with the reference template profile in 
the frequency domain (Taylor 1992). The time offset 
corresponding to each of the phase shifts found in the 
correlation was then added to the time-stamp recorded 
for each profile, resulting in a TOA that represents a time 
very close to the midpoint of each particular integration. 
In total, we measured 2446 TOAs: 106 from Parkes data, 
2233 from GBT data, and 107 from Nangay data. 

A model ephemeris was then fitted to the topocen- 
tric TOAs, using the tempo software package 3 . Included 
in this model is the motion of the Earth, calculated us- 
ing the JPL DE405 Solar System model (Standish 2004). 
To account for any instrumental and standard template 
profile differences, we fit for arbitrary time offsets for 
the Parkes and Nangay-derived TOAs, with reference to 
those from the GBT. For the GBT and Parkes data sets, 
corrections were also made to account for offsets between 
the clock readings from each observatory and UTC time, 
obtained using data from the Global Positioning System 
(GPS) satellites. In the case of Nangay data, recorded 
times are derived for UTC directly from GPS, and thus 
no additional clock corrections were needed. 

To obtain a best-fit value for dispersion measure (DM), 
we averaged the GBT-derived pulse profiles into five 
center-frequency bins (1348, 1364, 1384, 1404, and 1424 
MHz). We performed timing on this subset of the to- 
tal data set using the best-fit solution derived from all 
telescope data, allowing DM and its derivative to vary, 
while holding fixed all other system parameters. We ar- 
rived at a value for DM (149.6258 ± 0.0006 pecur 3 ) that 
we then held fixed for the timing analysis on the entire 
data set. We have also found evidence for the existence 
of a dispersion measure derivative, which we include in 
our timing model (see Table 2). In one observation (2006 
December 13), excess time delay was observed because 
the pulsar-Earth line-of-sight passed near the Sun, caus- 
ing a temporary increase in electron column density (e.g., 
Splaver et al. 2005; You et al. 2007). To account for this, 
we included an arbitrary time offset for this day as a 
parameter to be fit in the timing analysis. 

The effects of orbital motion on the pulse arrival times 
were taken into account using the ellI timing model 
(Lange et al. 2001), due to the near-circularity of the 
pulsar's orbit. This model parametrizes the eccentric- 

3 http : //www . atnf . csiro.au/research/pulsar/tempo 



4 



Fcrdman et al. 




0.4 o.e 
Orbital phase 

Fig. 2.— GBT-derived timing residuals for the PSR J1802-2124 
system, plotted against orbital phase relative to ascending node 
passage. For clarity, we have averaged the residuals into 256 orbital 
phase bins. Top: The full effect of Shapiro delay. Here, the Shapiro 
delay r and s parameters were excluded from the fit, with the best- 
fit orbital and other parameters held fixed. Middle: Once again, 
the Shapiro delay terms were left out of the fit, but in this case the 
Keplerian orbital parameters were left to vary as free parameters. 
Some of the Shapiro delay signal is absorbed into these parameters. 
However, the effect of Shapiro delay is still very evident in these 
residuals. Bottom: All parameters, including Shapiro delay, were 
included in the timing model fit. 



ity e and longitude of periastron uj in terms of the two 
parameters r\ = esinu; and k = ecoscj, which are used 
in the timing fit. It also replaces the time of periastron 
passage T with the time of ascending node T asc as the 
reference epoch to be fit in the timing model. 

In addition to the basic Keplerian parameters, we fit- 
ted for the Shapiro delay of the pulsed emission as it tra- 
versed the gravitational potential well of the companion 
star. This effect is described in the timing model in terms 
of the so-called "range" (r) and "shape" (s) parameters; 
the delay in the pulse arrival times for small-eccentricity 
orbits is given by: 



At 



-2r In < 1 — s sin 



2tt 
Pb 



(t - T a; 



(1) 



where t is the pulse TOA. Unless the orbit is close to 
edge-on, the Shapiro delay cannot be disentangled from 
the arrival time delay due to orbital motion (see Lange 
et al. 2001, appendix). Figure 2 shows the timing resid- 
uals over orbital phase, resulting from various fits to the 
TO As. The effect of Shapiro delay is still very evident 
when fitting for the Keplerian orbital parameters, which 
absorb some, but not all, of the Shapiro delay signal. 
Once the r and s parameters are measured, they can 
be converted into the companion mass mi and inclina- 
tion angle i. This is done through the following relations 
(Damour & Deruelle 1986): 

r=^ (2) 

s = sinz. (3) 

where the relation for r assumes that general relativity 
(GR) is the correct description of gravity, with m 2 ex- 
pressed in solar masses. 

Figure 3 shows the timing residuals from each instru- 
ment over time. In obtaining a best-fit model using these 
values, the scatter in the resulting residuals, while very 



close to having a random Gaussian distribution about 
zero, was larger than most of the errors on the indi- 
vidual data points, which were derived from the tem- 
plate profile cross-correlations. This resulted in an over- 
all value of x 2 per degree of freedom v that is greater 
than one (x 2 jv = 1.18 for Parkes and GBT data, and 
1.47 for Nangay). This was almost certainly due to an 
underestimation of the TOA uncertainties that resulted 
from the profile cross-correlation process, or from lower- 
quality profiles that arose because of signal contamina- 
tion by radio frequency interference, or perhaps coarse 
signal quantization as the data were sampled. To com- 
pensate, we have calculated an amount to add in quadra- 
ture to the original uncertainties in the TOAs, so that 
X 2 jv = 1.0 for each telescope data set. The GBT TOAs 
had errors that required very little correction, and dom- 
inated the data set. We thus, unless otherwise noted, 
report the uncertainties directly output by tempo as 68% 
confidence limits on the fit parameters, shown in Table 2. 
The individual weighted rms values of the post-fit tim- 
ing residuals for this pulsar are 6.1 /is from Parkes data, 
2.2 /xs from the GBT data, and 3.6 ^.s from the Nangay 
data. The combined value for the weighted rms of the 
residuals is 2.3 fis. Another compensation method that is 
used to arrive at x 2 l v = 1-0 involves the multiplication of 
the TOA uncertainties by calculated scaling factors. In 
doing this, we find a weighted rms of the post-fit residu- 
als of 2.2 /is, an approximately 5% lower value. However, 
we report the slightly more conservative parameter mea- 
surements found by using the former method; these can 
be found in Table 2. 

It should be noted that we obtain significant mea- 
surements of several spin frequency derivatives from 
our timing analysis (see Table 2). This may be at- 
tributed to intrinsic pulsar timing instability, though this 
is not typically seen in recycled pulsars — including those 
few with similar rotation periods and surface magnetic 
fields to PSR J1802-2124 (e.g., PSRs J0900-3144 and 
J1804-2717; see Burgay et al. (2006) and Hobbs et al. 
(2004), respectively) — with some exceptions (e.g., Kaspi 
et al. 1994; Matsakis et al. 1997; Verbiest et al. 2009). It 
may also be due to unmodclcd effects from the interstel- 
lar medium along the direction of PSR J1802— 2124, or 
from intervening material within the Solar System that is 
detectable as a result of the low ecliptic latitude of this 
pulsar. While fitting for these higher-order frequency 
derivatives (so-called "polynomial whitening"; see, e.g., 
Lorimcr & Kramer 2005) may affect the measurement of 
astromteric parameters such as position and proper mo- 
tion (reflected in Table 2 by quoting 2<r uncertainties), 
this long-term trend did not have a significant effect on 
our measurements of the orbital parameters, which is our 
focus in this work. 

4. RESULTS 

In order to ensure that the measured system masses 
represent the best model fit, we probed the x 2 over a fine 
grid of values, evenly distributed in |cosi|— mi space, and 
allowing all other timing parameters to vary. (We use the 
absolute value of cosi since we cannot distinguish i < 
90° from i > 90°.) The resulting confidence contours 
mapped by the x 2 values within our grid are shown in 
Figure 4. Overlaid are curves of constant pulsar mass mi, 
which is calculated through the Keplerian mass function, 
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TABLE 2 

Measured and derived parameters for PSR J1802— 2124 



Measured parameters 



Data span (MJD) 52605.2 - 54950.4 

Right Ascension, a (J2000) 18 h 02 m 05f 335576(9) (2a) 

Declination, 8 (J2000) -21°24'03'.' 649(3) (2a) 

Proper motion in a, fi a (mas yr _1 ) —0.85(10) (2a) 

Proper motion in <5, \fis\ (mas yr _1 ) < 4.8 (2a) 

Rotation frequency, v (s^ 1 ) 79.066424229950(2) 

Rotation frequency derivative, i> (10~ 16 s~ 2 ) —4.5360(16) 

Reference Epoch (MJD) 53453.0 

Dispersion measure, DM (peem -3 ) 149.6258(6) 

Derivative of dispersion measure, DM (10 — 5 peem -3 yr _1 ). . . — 9(3) 

Projected semimajor axis, x (lt-s) 3.7188533(5) 

Orbital period, P b (days) 0.698889243381(5) 

Epoch of ascending node passage, T asc (MJD) 53452.633290841(4) 

7; = esinu; 0.00000086(9) 

KEecosu 0.00000232(4) 

Cosine of inclination angle a , |cosi| 0.176(11) 

Companion mass a , mi (M Q ) 0.78(4) 



High-order rotation frequency derivatives 

Second frequency derivative (10 — 25 s~ 3 ) —1.04(7) 

Third frequency derivative (10 -32 s~ 4 ) —1.15(11) 

Fourth frequency derivative (10 -40 s~ 5 ) 9.2(5) 

Fifth frequency derivative (10~ 47 s~ 6 ) 2.4(4) 

Sixth frequency derivative (10 — 54 s~ 7 ) —4.7(4) 

Seventh frequency derivative (10 1 s~") 1.45(11) 



Derived parameters 

Rotation period, P (s) 0.0126475935865227(3) 

Rotation period derivative, P (lO^'ss^ 1 ) 7.256(3) 

Distance to pulsar b , d (kpc) 2.94 

Characteristic age, t c (Gyr) 2.76 

Surface magnetic field, B Bur f (10 8 G) 9.7 

Eccentricity, e (xl0~ 6 ) 2.48(5) 

Longitude of periastron, u) (°) 20(2) 

Epoch of periastron passage, T (MJD) 53452.673(4) 

Orbital inclination 21 , i (°) 79.9 or 100.1(6) 

Mass function, / (M ) 0.11305589(5) 

Pulsar mass a , m\ (M©) 1.24(11) 



Note. — Unless otherwise noted, parentheses indicate the la uncertainties on the last 
digit (or last two digits, if two digits are given). 

a Reported mass and |cos i\ values represent the median value from the respective probability 
density functions. 

b The distance is estimated using the NE2001 Galactic free electron distribution model 
(Cordes & Lazio 2002). 



given by: 
f(m 1 ,m 2 ,i) 



(m 2 sini)" 



2tt 



1 



(4) 



(mi+m 2 ) 2 \P b/ 

where x = asini/c is the projected semi- major axis of 
the pulsar's orbit, and T Q = G M Q / c 3 = 4.925490947 x 
1CP 6 s is the mass of the Sun in units of seconds. The 
mass function / and system masses mi and m 2 are ex- 
pressed in units of solar mass. 

The most probable pulsar mass, companion mass, 
and | cos i | were found by calculating their respective 
marginalized probability density functions (PDFs). De- 
tails of the method can be found in Splaver ct al. (2002, 
appendix). The intervals representing the 68% uncer- 
tainty in these quantities were calculated by determining 
the parameter values to either side of which the tails 
of each PDF cover 16% of the total area under the re- 
spective functions. We find that the best-fit median 



pulsar and companion masses are 1.24 ± 0.11 M Q and 
0.78 ± 0.04 Mq, respectively. This represents the second- 
most precise timing measurement to date of a pulsar in 
a NS-WD binary, after PSR J1909-3744 (Jacoby et al. 
2005). We note that the weighted rms of the post-fit 
timing residuals varies by less than 10% when the sys- 
tem masses and inclination angle are fixed at ±3cr from 
their best-fit values. 

We have also measured the right-ascension component 
of the system's proper motion to be fi a = —0.85 ± 
O.lOmasyr^ 1 (95% uncertainty). The pulsar's small 
ecliptic latitude has made it difficult to measure its 
proper motion in declination with the current data set. 
We thus quote a 95% upper limit of 4.8 mas yr" 1 for this 
quantity. However, we can calculate a one-dimensional 
space velocity based on our measurement of fi a and 
estimated distance of 2.94 kpc, based on the NE2001 
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Fig. 3.— Timing residuals for the PSR J1802-2124 system 
plotted as a function of time (Modified Julian Date). Parkcs, GBT, 
and Nancay data are represented by open squares, filled circles, and 
open diamonds, respectively. 



Galactic free electron distribution model (Cordes & Lazio 
2002), given the measured pulsar DM. We find that the 
pulsar velocity in the right ascension direction v a =12 
kms -1 , suggesting a relatively low velocity compared 
with those of other millisecond and binary pulsars (see, 
e.g., Toscano et al. 1999; Lommen et al. 2006). See §7 
for further discussion of the pulsar velocity. 

5. EVOLUTION OF THE PSR J1802-2124 SYSTEM 

The timing results from observations of 
PSR J1802— 2124 over the past six and a half years 
show that it is in a relatively compact binary system 
with a massive WD. It is also a light pulsar; along with 
others such as PSR J1141-6545 (Bailes et al. 2003) 
and PSR J1713+0747 (Splaver et al. 2005) for example, 
it is among the least massive known NSs with WD 
companions. This mass measurement represents the 
first made for what we refer to as short orbital-period 
IMBPs (P b < 3 days; see Table 3). 

We find that several NS-WD binary formation sce- 
narios cannot explain the observed parameters of the 
PSR J1802-2124 system. The usual LMBP mass trans- 
fer scenario, which would invoke an extended, stable pe- 
riod of accretion of matter onto the NS surface, is dif- 
ficult to reconcile with our measurements of the pul- 
sar and WD companion masses; it also seems to be in- 
compatible with the measured pulsar rotation period, 
which is significantly longer than those typically found in 
LMBP systems. The highly super-Eddington accretion 
scenario outlined earlier (Tauris et al. 2000) also does 
not appear to be able to produce the PSR J1802— 2124 
system — while the companion is likely to be a CO WD 
(0.4 < m 2 < 0.9 M Q ), the 16.8-hour orbital period is 
significantly less than the <~ 3-day minimum period pro- 
duced in this scenario. The double-He core progeni- 
tor scheme of B01, used to explain the formation of 
PSR B0655+64, assumes that the progenitors of the NS 
and WD have similar masses. However, the NS mass in 
the PSR J1802-2124 system differs significantly from 
that of the WD companion. This probably indicates 
a corresponding disparity in mass for their progenitors, 
making this theory difficult to apply in this case. This 
leads us to believe that the most viable formation sce- 
nario for the PSR J1802— 2124 system is that it had sur- 
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Fig. 4. — 68% (center), 95%, and 99.7% confidence contours 
for PSR J1802— 2124 in orbital inclination-companion mass space. 
Curves of constant pulsar mass are plotted over the contours. 



vived a phase of CE evolution. This is supported by the 
compactness of the orbit, the large WD mass, and mod- 
erately slow spin rate. Our mass measurements indicate 
that the pulsar probably had little time to accrete matter 
from the companion star before the envelope was ejected 
(~ 10 3 years; e.g., Tauris & van den Heuvel 2006), as 
its mass is similar to those of recycled pulsars in DNS 
systems (e.g., Nice et al. 2008; Stairs 2008) 



6. COMMON-ENVELOPE SURVIVAL RATE OF IMBPS 

Although we believe the evolution of the 
PSR J1802— 2124 system to be generally understood, 
the precise histories of this system and other short 
orbital-period IMBPs can only be elucidated through 
detailed binary stellar evolution simulations, which are 
well beyond the scope of this paper. We now briefly 
discuss in broad terms several potential contributions 
to the discordance between the theorized and observed 
population of IMBPs. 

6.1. Observational clues 

Using the population synthesis results of Porte- 
gies Zwart & Yungelson (1998), as well as evolution anal- 
ysis by Bethe & Brown (1998), B01 argue that, within a 
factor of two, the ratio of the birthrate of short orbital- 
period IMBPs (i.e., those which have passed through a 
CE phase in their evolution) to that for the young ec- 
centric NS-massive WD binaries is expected to be ~ 1. 
Only two of the latter class have been observed (Stokes 
et al. 1985; van Kerkwijk & Kulkarni 1999; Kaspi et al. 
2000a). B01 also argue that due to their longer visible 
lifetimes (due to their approximately two orders of mag- 
nitude weaker surface magnetic fields), we should observe 
a factor of ~ 50 — 100 more short orbital-period IMBPs, 
while only four or five have so far been seen. According to 
this argument, hypercritical accretion-induced collapse of 
most CE-embedded NSs into black holes is responsible 
for this observational discrepancy. 

We also note that through evolutionary model analysis, 
Belczynski et al. (2002) find that roughly two-thirds of 
NS systems can survive hypercritical accretion for a max- 
imum NS mass greater than 2M Q . However, we find a 
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relatively low mass for PSR J1802— 2124, indicating that 
this pulsar (as well as those in the precisely-measured 
pulsars in DNS systems) has not accreted a significant 
amount of matter, regardless of the specific formation 
mechanism undergone by these systems. 

The wide range of observed 1400-MHz luminosities for 
both short orbital-period IMBPs and eccentric NS-WD 
systems make a rigorous comparison of the expected and 
observed populations difficult without exploring an array 
of population synthesis models (e.g., Kim et al. 2004; 
Kalogera et al. 2005) or considering the various survey 
selection effects, some of which we now outline. 



6.2. Selection effects 

We present here two important selection effects against 
the observation of IMBPs that have been neglected by 
B01. Firstly, for a given luminosity and DM, a shorter 
spin period will render the pulsar more prone to the ob- 
servational effects of dispersive smearing. This is be- 
cause, in faster-rotating pulsars, the pulse will become 
smeared to a greater extent as a fraction of the spin 
period. This applies here, since IMBPs show an over- 
all spin period distribution that is substantially shorter 
than in the eccentric NS-WD binaries (see Table 3). In- 
deed, PSR J1802-2124 has a spin period-to-DM ratio 
that is approximately a factor of 40 smaller than, for 
example, PSR J1141-6545 (P spin /DM - 0.084 and 3.4, 
respectively). To further illustrate this point, we have 
calculated the observed fractional pulse width as a func- 
tion of DM for two pulsars at 1400 MHz, assuming 3- 
MHz channels. This is similar to the search observa- 
tion setup for the Parkes Multibeam Pulsar Survey (e.g., 
Manchester et al. 2001). An eccentric NS-WD binary 
pulsar like PSR J1141-6545 (P spin = 394 ms; Kaspi et al. 
2000a) has a fractional pulse width that would never 
be smeared by more than 2% of the pulse period, out 
to a DM > 1200 pccm~ 3 , near the limit of the known 
pulsar population; hence dispersion smearing is negligi- 
ble for this class of pulsars. In contrast, an IMBP like 
PSR J1802-2124 (P spin = 12.6 ms) has a observed frac- 
tional pulse width that grows to more than one-sixth 
of a pulse period, making the pulsar relatively difficult 
to detect, at DM ~ 190pccm~ 3 or larger. This im- 
plies that the eccentric young-pulsar binaries can be dis- 
covered to much larger volumes than the IMBPs. The 
largest DM to which an IMBP has thus far been discov- 
ered is PSR J1810-2005, at 240.2 pccuT 3 (Camilo et al. 
2001), which has a spin period approximately three times 
that of PSR J1802-2124. 

To quantify the relative detectability and survey vol- 
ume due to propagation effects, we performed the fol- 
lowing simple Monte Carlo population using the freely 
available psrpop pulsar population modelling package 
(Lorimer et al. 2006a). We generated model galaxies 
containing pulsars with periods and pulse widths iden- 
tical to PSR J1802-2124 and PSR J1141-6545. Each 
model population is distributed in Galactocentric radius 
according to the best-fit distribution found by Lorimer 
et al. (2006a). To test whether there is any dependence 
on the dispersion above the Galactic plane, we gener- 
ated models with exponential z-height distributions with 
means in the range 300-500 pc. Since we are only inter- 



ested in propagation effects on the relative detections of 
the two pulsars, we assign each model pulsar a radio lu- 
minosity of 100 mJy kpc 2 . For each population, we then 
record the number of detected pulsars in the detailed 
model of the Parkes Multibeam Pulsar Survey described 
by Lorimer et al. (2006a). Regardless of the assumed z- 
scale height, the ratio of detections of PSR J1141— 6545 
to PSR J1802-2124 is always ~ 2.2. In other words, 
due to propagation effects alone, the effective volume of 
the Galaxy surveyed by the Parkes Multibeam Survey 
for PSR J1802-2124-like objects is only half that of pul- 
sars similar to PSR J1141— 6545. It should be noted 
that this simulation does not take account of any beam- 
ing corrections which are likely to be significant given 
that longer period pulsars generally have smaller beams 
than their shorter period counterparts (e.g., Tauris & 
Manchester 1998). To quantify this, using the Tauris & 
Manchester (1998) beaming model, we estimate that the 
beaming fraction of PSR J1802— 2124 is four times larger 
than that of PSR J1141-6545. When combined with the 
above simulations, this would imply that for populations 
of comparable sizes, one might expect only half as many 
objects that are PSR J1141— 6545-like compared to those 
that are similar to PSR J1802— 2124. Future surveys that 
use narrower-channel instruments should expect to find 
a larger number of higher-DM IMBPs compared to ec- 
centric NS-WD binaries. 

We also note that pulsar acceleration would cause fur- 
ther signal spread in the Fourier search domain, to a 
larger extent for IMBPs than for young pulsars in sim- 
ilar orbits (e.g., Hessels et al. 2007). Although difficult 
to quantify, this is an important additional selection ef- 
fect against discovery of IMBPs relative to the eccentric 
NS-WD binaries. 

While the discrepancy in the observed ratio of these 
two system types pointed out by B01 may still be sup- 
ported by the available data on IMBPs, we emphasize 
again that a more precise estimation of observable num- 
bers of short orbital-period IMBPs will only come with 
further population synthesis studies, as well as account- 
ing carefully for survey selection effects. 

7. FUTURE MEASUREMENTS AND STUDIES 

A more precise measurement of the full proper mo- 
tion will help us to constrain the space velocity of this 
system. This will further our understanding of IMBP 
formation history. For example, it would allow us to test 
the apparently low scale height of IMBPs compared to 
LMBPs, presumably due to the larger combined mass of 
IMBP progenitor systems, as suggested by Camilo et al. 
(2001). The low transverse velocity that is hinted at by 
the measured proper motion in the right ascension direc- 
tion supports this theory. Within five years, we expect 
to obtain a significant measurement of proper motion in 
declination. This will enable us to calculate a more reli- 
able value for the transverse velocity of PSR J1802— 2124. 

Based on the measured orbital parameters, the GR 
prediction for orbital decay P & GR = —3.1 x 10~ 14 ss _1 
for the PSR J1802— 2124 system, and we expect to sig- 
nificantly measure p, in the near future. This has only 
been achieved in four NS-WD systems: in the case of 
PSR J1141-6545 (Kaspi et al. 2000b; Bailes et al. 2003; 
Bhat et al. 2008) and PSR J0751+1807 (Nice et al. 2005, 
2008), the p, intrinsic to the system is measured; in the 
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TABLE 3 

Known short orbital-period intermediate-mass binary pulsars and eccentric neutron star-white dwarf systems 



PSR 


Rotation 


Orbital 


Eccentricity 


Surface Flux density 


Distance 3, 


Luminosity density 


References 




period 


period 




magnetic field (1400 MHz) 


(kpc) 


(1400 MHz) 






(ms) 


(days) 




(10 10 G) ' (mjy) 




(mjy kpc 2 ) 




Short orbital-period IMBPs (P b < 3 days) 


B0655+64 


195.7 


1.03 


0.0000075 


1.17 0.3 


0.49 


0.9 


1,2 


J1232-6501 


88.2 


1.86 


0.00011 


0.856 0.34 


6.2 


160 


3,4 


J1435-6100 


9.34 


1.35 


0.0000105 


0.0484 0.25 


2.2 


15 


3,4 


J1757-5322 


8.87 


0.453 


0.0000040 


0.0489 2.3 


0.96 


27 


5,6 


J1802-2124 


12.6 


0.699 


0.0000025 


0.0966 0.77 


2.9 


84 


7, this work 


Eccentric NS-CO WD binaries 


J1141-6545 


394 


0.198 


0.17 


132 3.3 


> 3.7 b 


> 570 


8,9 


B2303+46 


1066 


12 


0.66 


78.8 0.38 c 


2.9740 


42 


10 



References. — 1. — Jones & Lyne (1988); 2. — Lorimer et al. (1995); 3. — Camilo et al. (2001); 4. — Manchester et al. (2001); 
5. — Edwards & Bailes (2001b); 6 — Bailes, private communication; 7. — Faulkner et al. (2004); 8. — Kaspi et al. (2000a); 9. — Ord 
ct al. (2002); 10.— Stokes et al. (1985). 



a For all pulsars, except where noted, we have used the NE2001 model (Cordes & Lazio 2002) to derive the distances, based on 
the dispersion measures of the objects. 

b For PSR J1141— 6545 we have used the distance published in Ord et al. (2002), derived from the neutral hydrogen absorption 
spectrum along the pulsar's line of sight. 

C PSR B2303+46 does not have a published 1400 MHz flux density. In this case we estimate it using published spectral index for 
this pulsar, from Maron ct al. (2000). 



PSR J1012+5307 system (Lazaridis et al. 2009), the in- 
trinsic Pf, is not yet significantly determined, however 
the contribution to the measured value from kinematic 
effects is of the same order as the quadrupolar GR pre- 
diction; and in PSR J0437-4715 (Vcrbicst et al. 2008), 

the measured Pb is attributed predominantly to contam- 
ination by apparent acceleration due to the system space 
velocity, as the measured Pf, is four orders of magnitude 
greater than the value predicted by GR for that system. 

For PSR J1802— 2124, we anticipate an intrinsic i\ 
measurement which, when combined with the Shapiro 
delay measurements, will overdetermine the system mass 
values and allow for a check on the relativistic analysis 
of the kinematics of the system. 

The difference in self-gravities between 
PSR J1802-2124 and its WD companion is pre- 
dicted to cause a deviation in Pb from the quadrupolar 
GR prediction according to some scalar-tensor theories 
of gravity (e.g., Esposito-Farese 2005). A limit to this 
departure from the GR value could be used to constrain 
the existence of dipolar radiation from this system, or to 
set a limit to the variation of the gravitational constant, 
G (Nordtvedt 1990; Damour & Taylor 1991; Lazaridis 
et al. 2009). It remains to be seen, however, to what 
extent the kinematic corrections (Damour & Taylor 
1991; Nice & Taylor 1995) may affect our ability to use 
this measurement for any of these purposes. 

8. CONCLUSIONS 

In this paper we have described the observations and 
timing analysis of PSR J1802-2124. These have pro- 
vided updated system parameters and precise measure- 
ments of the pulsar and companion WD masses, which we 
find to be 1.24±0.11 M and O.78±O.O4M , respectively. 
In particular, this determination was made possible by 
the detection of Shapiro delay on the pulse arrival times. 



The result is of particular significance for this class of pul- 
sar, since it is only the second such mass measurement 
for an IMBP system, and the first for a short orbital- 
period IMBP (Pb < 3 days). The mass measurements of 
the PSR J1802-2124 system highlight the dependence 
of the final system configuration on the specific mass- 
transfer history and particular evolution of the system 
in question. The similarity between the properties of 
this system to those of recycled DNS pulsars hints that 
the evolutionary paths of these two system types may be 
analogous, probably involving a CE/inspiral phase. 

It is clear that to arrive at a definitive picture of 
IMBP evolution, and more generally, the evolution of 
the many observed binary system types, we must dis- 
cover more systems with measurable masses, and several 
recent and future surveys (e.g., van Leeuwen & Stappers 
2004; Cordes et al. 2006; Lorimer et al. 2006b) are ex- 
pected to find many such systems for study. 
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